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Abstract 

We describe some numerical experiments which determine the degree of 
spectral instability of medium size randomly generated matrices which are 
far from self-adjoint. The conclusion is that the eigenvalues are likely to 
be intrinsically uncomputable for similar matrices of a larger size. We also 
describe a stochastic family of bounded operators in infinite dimensions for 
almost all of which the eigenvectors generate a dense linear subspace, but 
the eigenvalues do not determine the spectrum. Our results imply that the 
spectrum of the non-self-adjoint Anderson model changes suddenly as one 
passes to the infinite volume limit. 
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1 Introduction 

In a series of recent papers [2,4-10] a number of authors have investigated the spec- 
tral properties of non- self-adjoint matrices and operators, coming to the conclusion 
that the eigenvalues are frequently highly unstable under small perturbations of 
the coefficients of the matrices. Trefethen has investigated a series of numerical 
examples using the concept of pseudospectrum (the contour plots of the resol- 
vent norm), which provides a graphical demonstration of the degree of instability, 



20, 21]. This concept is, however, not well adapted to the consideration of very 



large numbers of randomly generated matrices, for which one needs to produce a 
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numerical measure of instability of the spectrum. In this paper we define such an 
instability index, and compute it for a series of randomly generated N xN matrices 
for various values of iV up to 50. Our numerical results are presented in Section 
3, following a short theoretical section which describes the concepts involved. Our 
results are fully in line with what would be expected by experts in pseudospectral 
theory, but we believe that such a systematic quantitative investigation of non-self- 
adjoint random matrices has not previously been carried out, and it is clear that 
much of the spectral theory community is not aware of these phenomena. 

From Section 5 onwards we consider a related problem for a stochastic family of 
non-self-adjoint bounded operators in infinite volume, i.e acting on Z 2 (Z n ). For 
one example we prove that although the eigenvectors of almost all of the operators 
span a dense linear subspace, they do not form a basis, and the spectrum is much 
larger than the closure of the set of eigenvalues. Section 6 is devoted to spelling out 
the implications of our results for the non-self-adjoint Anderson model of Hatano- 
Nelson, which has been the focus of much recent attention. [12, [13] , [14], |17| We find 
that the asymptotic behaviour of the eigenvalues as the volume increases does not 
describe the full spectrum of the infinite volume problem. The reason is that there 
are many approximate eigenvalues of the finite volume problem which are not close 
to true eigenvalues but which nevertheless affect the infinite volume limit. This 
situation is typical of non- self-adjoint operators, for which the eigenvectors need 
not be even approximately orthogonal. 



2 The Theoretical Context 

Throughout this section we suppose that A is an iV x N matrix with distinct eigen- 
values. This is generically true, since the set of matrices with repeated eigenvalues 
forms a lower dimensional set with zero Lebesgue measure. All vectors below are 
assumed to be column vectors unless otherwise stated, and * denotes the conju- 
gate transpose. Norms of vectors are always their Euclidean norms and norms of 
operators are always the corresponding operator norms. 

Let A be an eigenvalue of A of multiplicity 1, with corresponding eigenvector <ft. Let 
if) be an eigenvector of A* with eigenvalue fi. If fi ^ A then ip*<f) = and otherwise 
this inner product is non-zero. Assuming that = A, the spectral projection P 
associated to the eigenvalue is given by 

pf = {r^irm 

and its norm is 

\\p\\ = uwwmmv w 

If we label all the above quantities by n for n = 1, N, then we may define the 
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instability index of A by 

i(A) = max{||P n || :n=l, N}. 
This is unrelated to its condition number 

k(A) = \\A\\ WA^W 

since the instability index of a self-adjoint or normal matrix always equals 1, but 
its condition number may be arbitrarily large; in the converse direction the matrix 




has small condition number but unbounded instability index as 8 — > 0. The norm 
of a particular spectral projection is also called the condition number of the eigen- 
value, and is known to measure how unstable the eigenvalue is under small pertur- 
bations of the matrix, [[l], |23|. 24] . [16, sect. 11.2]. We emphasize that if the norm 



of the spectral projection is very large the instability of the eigenvalue it intrinsic: 
it does not depend on the particular method of computing it. The norm of P is 
always at least 1 and equals 1 if and only if P is orthogonal, or equivalently if 
(j) = if). The following proposition relates this index to other measures of how far 
the matrix is from being normal, Jl|, 23, 24 . 



Proposition 1 We have the relations 

(i) =>- (ii) =>- (Hi) =>- (iv) 
between the following conditions, the constant k being the same in all cases. 

(i) There exists an invertible matrix V with k(V) < k such that D = V~ 1 AV is 
diagonal. 

(ii) The functional calculus satisfies 

11/04)11 <*||/IU 
for all complex-valued functions f defined on Spec(A), where 

||/|| 0O = max{|/(A)|:AeSpec(A)}. 
(Hi) The resolvent operators satisfy 

|| (A - zl^W < k dist(z, Spec(A)) -1 

for all z ^ Spec(v4) ; where dist is the Euclidean distance of a point from a set. 

(iv) The spectral projection P of every eigenvalue X of A satisfies \\P\\ < k, and 
hence 

i(A) < k. 
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Proof 

(i) (ii) We use f(A) = Vf(D)V~\ Spec(A) = Spec(L>) and 

11/(0)11 = ll/lloo. 

(ii) =>■ (iii) This is a matter of considering the particular function /(A) = (A — z)~ l . 

(iii) =>- (iv) We express the spectral projection as a contour integral of resolvent 
operators around a small circle centred at A. 

Note 2 The proof of the theorem remains valid if we replace the use of the Eu- 
clidean norm on by any other norm, provided the appropriate operator norm 
is used for matrices and the operator norm of any diagonal matrix D is given by 

||D|| = max{\D n>n \ : 1 < n < N}. 

This is equivalent to the norm being absolute sect. 10.5] and holds in particular 
for all of the V norms. However, certain other matters, such as the identification 
of orthogonal projections with those of norm 1, are dependent on the use of the 
Euclidean norm. 

Note 3 If we assume condition (iv) of the above theorem then it follows from the 
formula 

(a — z)~ l — Y^Pi(K — zy 1 

i 

that 

II (A - z)- 1 1| < Nk dist(Spec(A), z)' 1 
Since the value of k frequently increases exponentially with the dimension iV,pO" 



p25f one must expect pseudospectral information and that obtained from the insta- 
bility index to be broadly equivalent. 

The point of the theorem is that if any spectral projection of A has very large 
norm, then the constant k of any of the earlier conditions must be very large, and 
diagonalization of the matrix A is an intrinsically ill-conditioned procedure. 

There is a family of matrices for which the instability index defined above can be 
computed in closed form. This is of some interest for its own sake, but we used it to 
verify the algorithm used to compute the instability indices of randomly generated 
matrices. We assume that 

A m , n = f(m-n)a m - n 
where a > 1 and / : Z — > C is any function which is periodic with period N. 
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Theorem 4 The eigenvectors of A are of the form 

<P r (n) = a n e 27rirn/N 



where r = 1, N. The corresponding spectral projections all have norm 



(a 2 



Thus the instability index of A is also c. 

Proof The first statement is a matter of applying the matrix to such a vector, 
and noting that the set of all such vectors is a basis for C N . The corresponding 
eigenvectors of A* are 

^ r ( n ) = a - n e 2nirn/N . 

From these facts it is now easy to calculate the condition numbers of each of the 
eigenvalues using ([[]). 



3 Numerical Results 

We have applied the above ideas to a series of randomly generated tridiagonal 
N x N matrices for various values of N. The matrices are of the form A mj1l so 
that A m>n = for \m — n\ > 1, the other coefficients being chosen randomly and 
independently. If m — n = 1 the coefficients are chosen using a uniform distribution 
on [0, 1], if m — n = the coefficients are chosen using a uniform distribution on 
[0,2], and if m — n = — 1 the coefficients are chosen using a uniform distribution 
on [0,3]. 

For each A we used Matlab to obtain an invertible matrix V and a diagonal matrix 
D such that V~ 1 AV = D. The columns of V are then the eigenvectors <p n of A, 
and are provided by Matlab in normalized form. The rows of V^ 1 are i\)* m where 
ip n are the eigenvectors of A*. We get = 1 automatically, so 

lip || 2 = \\ib || 2 = E 

II 1 n || || Vn || - t -'n,n 

where E = V" 1 (V^ 1 )* . We thus finally get 

i(A) = max{|£ njn | 1/2 : 1 < n < N}. (2) 

For each of M randomly generated N x N matrices A we computed the instability 
index using (0), and then sorted the data points into increasing order. Defining 
V r to be the number such that such that r% of the instability indices were less 
than V r , we determined V50 and P95 for various values of N. In fact we carried 
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out each computation twice in order to give some idea of the degree of reliability 
of our results; we tabulated the average of the two values and the difference V r 
expressed as a proportion of the average. 



N 


M 


^50 


^50 


^95 


^95 


log(P 50 )/iV 


log(P 95 ), 


10 


10 6 


18.47 


0.007 


304.1 


0.001 


0.292 


0.572 


20 


10 6 


684.9 


0.005 


6.13 x 10 4 


0.008 


0.326 


0.551 


30 


4 x 10 5 


2.138 x 10 4 


0.013 


7.48 x 10 6 


0.001 


0.332 


0.528 


40 


10 5 


5.77 x 10 5 


0.047 


6.68 x 10 8 


0.050 


0.332 


0.508 


50 


10 5 


1.43 x 10 7 


0.016 


4.60 x 10 10 


0.032 


0.330 


0.491 



It is clear from our results that both P50 and P95 increase extremely rapidly with 
N. In fact our results support the conjecture that V50 ~ e N ^ 3 as N —> 00. Matlab 
is already having a little difficulty in computing the eigenvalues of the matrices for 
N = 50, and fails entirely for iV = 200. Such exponential increase has also been 
found by Trefethen et al for other models using pseudospectral methods. [25[ 

An alternative approach to the above questions would be to compute the expected 
value of i(A) over a large sample of matrices A, but this would have the disadvan- 
tage of being unduly influenced by the very large size of i(A) for a small proportion 
of choices of A. 



4 Distribution of Norms of Spectral Projections 

Instead of studying i(A) for randomly distributed matrices A, one may examine 
how the norms of the individual spectral projections are distributed. This may be 
done in several ways. In the first we sort the norms of the spectral projections of a 
particular N xN matrix A in increasing order, but instead of examining the largest 
of these, namely i(A), we evaluate the number j(A) half way through the list. We 
did this for a series of 10 5 randomly generated 30 x 30 matrices, and discovered 
that for 50% of these matrices j(A) < 221.2. The fact that this number is so much 
smaller than i(A) under the corresponding conditions indicates that only a small 
proportion of the spectral projections of a typical random matrix A have really 
large norms. 

A second procedure is to consider all of the norms of the spectral projections of 
the 10 5 randomly generated 30 x 30 matrices as one list of 30 x 10 5 numbers. 
When we carried out this numerical experiment we found that 50% of the norms 
so obtained were less than 157.8. There is no reason why the above two numbers 
should coincide, but they are of the same order of magnitude, leading to the same 
conclusion. 
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The final and most interesting method is to carry out a spectral analysis of the co- 
variance matrix associated with the norms of the spectral projections. We proceed 
as follows. 

For each matrix A we define the numbers X n > for 1 < n < N by 

X„ = log(||P n ||) (3) 

where these are sorted in increasing order. (It is possible to carry out similar cal- 
culations without taking the logarithm above, but the results are less compelling.) 
As A varies within the usual class these provide a family of N non-negative random 
variables whose covariance matrix is defined by 

The eigenvalues and eigenvectors of this matrix provide information about the 
distribution of the norms of the spectral projections of 'typical' random matrices. 

We carried out the above computation for a sample of 10 5 randomly distributed 10 x 
10 matrices. The eigenvalues of G were found to be 0.0060, 0.0196, 0.0239, 0.0323, 
0.0433, 0.0613, 0.1034, 0.2133, 0.7169, 54.6547. The fact that one eigenvalue is 
so dominant is very striking, and indicates that to a very good approximation 
most of the random matrices have very similar distributions of the norms of their 
spectral projections. We repeated the computation for a sample of 10 5 randomly 
distributed 30 x 30 matrices. There were only 4 eigenvalues larger than 1, these 
being 1.5, 3.1, 12.02, 1322.9. 

In both cases we computed the eigenvector v corresponding to the largest eigenvalue 
of the covariance matrix C. We found that v n was close to being proportional to 
n. This corresponds to the norms ||-P n || of a 'typical' random matrix A increasing 
exponentially with n. We conjecture that the dominance of the leading eigenvalue 
increases and the rate of increase of the spectral norms becomes more accurately 
exponential as the size N of the matrices increases. 

Let us be more precise about this. Given N we consider the sample space 

Q N = [O^f- 1 x [0,2]" x [0,3]^ 

provided with the uniform probability distribution. For each w £ fijv we described 
how to construct a matrix A^ and then the N random variables X n defined by (0) . 
These have a symmetric N x N covariance matrix C whose eigenvalues may be 
ordered so as to satisfy 

< Ai < ... < Atv- 

We finally define 

Ajv-i 
= -r — 
an 
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making explicit the dependence of fj, on N. The conjecture is then that 

lim /i N — 0. 

We tested this hypothesis by considering a series of T randomly generated N x N 
matrices for various values of T and N. The results below provide some support 
for the conjecture. 



N 


T 




Xn 


Hn 


10 


10 6 


0.725 


54.99 


0.0132 


20 


10 6 


4.422 


422.40 


0.0105 


30 


4 x 10 5 


12.02 


1319.5 


0.0091 


40 


10 5 


24.15 


2928.4 


0.0082 


50 


10 5 


40.73 


5408.9 


0.0075 



5 Operators with Randomly Distributed Coeffi- 
cients 

In this section we present some ideas relating to a random family of bounded linear 
operators acting on the infinite-dimensional Hilbert space / 2 (Z). This may be re- 
garded as the infinite volume limit of our earlier problems, although pseudospectral 
theory suggests that one should also study the 'same' operator on Z 2 (Z + ). Phys- 
ically the choice between these two operators depends upon whether one wishes 
to include end effects, which are present both for large finite intervals and for the 
operator on Z 2 (Z + ). Our methods can easily be adapted to this case, but we do 
not spell out the modifications needed. 

We first formulate the ideas at a moderately general level, and only later restrict 
attention to the non-self-adjoint Anderson model of Hatano-Nelson. Let A be the 
operator associated with an infinite matrix {A, m , n } where A mtH — if \m — n\ > 
1; we suppose that the vectors v n = (A n+ i >n , A n>n , A n>n+ i) e C 3 are distributed 
independently according to a common law //, where \x is a probability measure on 
C 3 with compact support K. It follows from the assumptions and the fact that 
the matrix {A mtH } is tridiagonal that it is associated with a bounded operator A 
such that 

3 

\\A\\ < max{|uy| : w G K}. 

r=l 

The above procedure defines a stochastic family of bounded operators A^, for u 
in the sample space K°°. It is of some interest that the results which we obtain 
are not truly probabilistic: the statements of our theorems only involve the set K 
rather than the probability measure \i. 
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In order to prove some results about the spectra of operators in the family we 
introduce a notion from [fTI |. Given any bounded operator X on a Hilbert space 
7i, we say that the operator Y lies in its limit class, Y G C(X), if there exists a 
sequence U s of unitary operators on TC such that U*XU S converges strongly to Y 
as s — > oo. We also define the approximate point spectrum o~(X) of X to be the set 
of all A G C for which there exists a sequence of vectors f s G TC such that ||/ s || = 1 
and 

lim \\Xf B -Xf a \\ = 0. 

In all the examples in this paper cr(X) = Spec(X) almost surely but it seems 
desirable for the sake of possible future applications to keep the logical distinction. 
The following two lemmas are modifications of a classical result of Pastur stating 
that the spectrum of random tridiagonal operators is almost surely constant; the 
proof uses translation ergodicity of the class of operators. @, p 167] 

Lemma 5IfY lies in the limit class of X then o~{X) I) °"00- Hence 

Spec(X) D \J{a(Y) : Y G C(X)}. 

In particular if each lies in the limit class of the other then a(X) = cr(Y). 

Proof Suppose that for some A G C and all e > there exists / G TC such that 
= 1 and \\Yf — A/|| < e. Now put f s = U s f and observe that 



\\Xf s - A/ s || = \\U:XUJ - Xf\\ -> \\Yf - \f\\ 

as s — > oo. Therefore \\Xf s — Xf s \\ < e for all large enough s, and A G o~(X). 
We apply the above to the randomly generated operators acting on Z 2 (Z). 

Lemma 6 The limit class of almost every operator A w generated as described con- 
tains every operator such that uj G K°°. 

Proof Let N G Z + and Cu = {v n } neZ G K°°. Given £>0we put 

V n = {z G C 3 : \\v n - z\\ < e} 

so that fJ,(V n ) > for all n. Taking U s : / 2 (Z) — > l 2 (Z) to be the unitary operators 
associated with appropriate translations of Z, we need to show that for almost 
every uj = {v n } G K°° there exists M such that v n+ M G V n for all —N < n < N. 

To prove this we put V = I\n=-N Vn and 

U = { v n}neZ = { w m}meZ 

where 

W m = {v m (2N+l)+r}r=-N e R 2N+1 . 
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The vectors w m are independent and identically distributed with positive proba- 
bility that w m G V. Hence the probability that none of the w m lie in V is zero. 

As an application of the lemma, let v , w G C 3 and let B V)W be the bounded operator 
associated with the infinite matrix {B m n \ such that B mn — if |m — n\ > 1; we 
also assume that the vector (B n+ i jTl , B Utn+ i) equals v if n > 1 and equals w if 
n < 0. The point of introducing this class of operators is that they are similar to 
operators whose spectrum is well understood. [0, [| |T5| ] 

Corollary 7 probability 1 one /ias 

Spec(A,) D (J {a(B %w ) :v,weK}. 

Proof Every operator B V)W lies in the limit class of A by Lemma 5. 

The significance of the above lemmas is best seen by applying them to an example. 
We assume that /i is a probability measure concentrated on a finite subset F of 
R 2 (more general probability measures can also be treated). We assume that 
/i(0, 0) > and that any other point (x, y) with fi(x,y) > satisfies x > and 
y > 0. We then define the random family of operators A on Z 2 (Z) as described 
above with the following simplification: we assume that all A n ^ n = and that the 
vectors A n are distributed independently according to the law of \i. 

Theorem 8 Depending to the choice of fi, either A^ is self-adjoint with probability 
1 or it is non-self-adjoint with probability 1. In the latter case A^ possesses a 
countable set of eigenvalues whose corresponding eigenvectors span a dense linear 
subspace of 'I 2 : (Z). With probability 1 the set of eigenvectors is not a basis o/7 2 (Z). 

Proof The self-adjoint case occurs when the support of \x is contained in the diag- 
onal set {(x, y) G R 2 : x = y} and we assume that this is not the case below. 

With probability 1 an infinite number of the pairs are equal to 

(0,0), and we assume that this happens for the increasing sequence {N r } where 
r G Z. It may then be seen that A w can be decomposed as the orthogonal direct 
sum of matrices C r of sizes M r x M r where M r = N r — iV r _i. Since each matrix 
C r is tridiagonal with positive off-diagonal entries and zero diagonal entries, its 
eigenvalues are all real and of multiplicity one. By combining all of the eigenvectors 
of the C r as r increases we see that the eigenvectors of A u almost surely span a 
dense linear subspace of Z 2 (Z). It remains only to prove that this set of eigenvectors 
is almost surely not a basis. 

Let (x,y) be any point in R 2 with < x < y and fi(x,y) > 0. Since the pairs 
of coefficients of A^ are chosen independently, among the C r there almost surely 
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exist all s x s matrices of the form B s where 

{x if m = n + l 
y if m = n — 1 
otherwise. 

Now the eigenvectors of each B s are given explicitly by 

TTi) \y) 

and the corresponding eigenvectors of the adjoint operator are 

where k labels which eigenvector is being considered and r which coefficient. 
It follows that 

\\M 2 

\m\ 2 

An application of (|IJ) now shows that lim^oo i(B s ) = +oo, from which it follows 
that with probability 1 the norms of the spectral projections of are not uniformly 
bounded. This implies that the eigenvectors cannot form a basis. 

An example of Zabzyk which has some similarities to those of the above theorem is 
discussed in Theorem 2.17 of 0, where the failure of the spectral mapping theorem 
for semigroups is demonstrated. Another type of example involving differential 
operators whose eigenvectors do not form a basis was presented in || ^ [TTJ. It 
appears that such a situation is relatively common for non-self-adjoint operators 
in infinite dimensions. Since the eigenvectors do not form a basis for the Hilbert 
space, there is no reason to expect that they determine the spectral behaviour of 
A^. We continue with the hypotheses formulated before Theorem 7. 

Theorem 9 The eigenvalues of the operator A u are almost surely all real. How- 
ever the spectrum of A^ almost surely contains the interior of the ellipse E defined 
by 

E = {xe ie + ye~ ie : 9 E [0,2tt]} 
for every non-zero (x,y) G F. 



J2 sin 

r=l 
s 

r=l 

s 

J2 sin 

r=l 



7T 



kr 



+ 
nkr 

Ttkr 



1 
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Proof The operator A w is almost surely the orthogonal direct sum of tridiagonal 
matrices C r , where we continue to use our previous notation. Each C r is similar 
to a real symmetric matrix and therefore has real eigenvalues. 

If v = (x, y) is a non-zero point in F then by Corollary 6 the spectrum of A u almost 
surely contains a(B 0jV ). We assume that < x < y, the case < y < x having a 
similar analysis involving the adjoint operators. A direct computation shows that 
if z lies inside the said ellipse then both solutions Wi of 

x — zw + yw 2 = 

satisfy \wi\ < 1. It follows that z is an eigenvalue of B 0iV , the eigenvector being 
the sequence / G / 2 (Z) given by 

f _ / w i ~ w 2 if r — 1 
r I otherwise. 



6 The Non-Self- Adjoint Anderson Model 

In this section we consider non-self-adjoint Anderson-type operators of the form 

Hf n = e"Vn-i + e 9 /n + i + V n f n (4) 

where g > and V is a random real-valued potential. We assume that the values 
of V at different points are independent and identically distributed according to a 
probability law \i which has compact support M C R. 

A considerable amount of attention has already been paid to such operators, which 



arise in population biology and solid state physics, as cited in [|12| [L3], [14], [L7j. If 
one supposes that the operator acts on l 2 {— N, N} subject to Dirichlet boundary 
conditions then its spectrum is almost surely real, because H is similar to the 
operator defined by the same formula but with g = 0. The similarity is determined 
by the operator 

Sf n = e n9 f n 

which is bounded on this finite-dimensional space. On the other hand if one imposes 
periodic boundary conditions, then the spectrum of H is much more interesting 
and has been analyzed in great detail in the limit of large iV both numerically and 
theoretically. 0, [13], [14], [17|] The imposition of periodic boundary conditions was 
justified in Nelson-Shnerb by the fact that they were considering a rotationally 
invariant problem (the diffusion of bacteria in a rotating nutrient). Goldsheid has 
observed that the spectral results proved in JT2|| for periodic boundary conditions 



are equally valid for a wide range of quasi-periodic boundary conditions. However, 
Trefethen has pointed out that the spectral behaviour of this type of operator 
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is highly problematical. |22| This phenomenon has been investigated from several 
points of view over the last decade, 0, || 0, [19|, |2(| ^TJ and among the conclusions 



is the warning that one cannot assume that a solution of a non-linear equation is 
stable simply because the eigenvalues of its linearization about the solution have 
negative real parts. 

We show that if one considers the non-self-adjoint Anderson model H acting on 
/ 2 (Z), the spectrum is entirely different from what one obtains by letting N — > oo 
subject to any of the above types of boundary condition. The following theorems 
are directed towards locating the spectrum of H, but it would clearly be highly 
desirable to find a precise formula for it. 

Theorem 10 Let H be defined on Z 2 (Z) by (fjj where V satisfies the stated condi- 
tions. Then 

Spec(H) C conv(E) + conv(M) 
almost surely, where conv denotes convex hull and E is the ellipse 

E = {e 9+ie + e~ 9 - ie : 9 E [0, 2tt]}. (5) 

Proof Let A be the operator obtained from H by deleting the potential V, so that 
Spec(A) = E. Then the result follows from the equation 

Spec(F) C Num(F) C Num(A) + Num(\/) = conv(£) + conv(M). 

where Num stands for the numerical range 

Num(A) = {(Af,f) : \\f\\ — 1}. 

Theorem 11 Under the above assumptions we also have 

Spec(iT) C{zeC: dist(z,£T) < m} 
almost surely, where dist is the distance function and 

m = max{|s| : s G M}. 

Proof Since A is a normal operator with spectrum equal to E, we have 

\\(z - Ay 1 \\ =dist(z,£)- 1 

for all z ^ E. Also || V]| = m almost surely, and the result follows by examining 
the convergence of the perturbation series 



oo 

r.-A-VY^Y.iz-A)- 1 (V{z-A)^) n . 

n=0 



In the reverse direction we may apply Lemma 4 and Theorem 5 to obtain the 
following result, which is further improved in Theorem ITB. 
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Theorem 12 With probability one, the spectrum of H contains the set E + m for 
every m G M. 



Proof The operator A + ml lies in the limit class of H for all m G M and its 
approximate point spectrum is E + m (as is its spectrum). 

Example 13 Let fi, the measure determining V, have support [—B,B]. By ap- 
plying Theorems ITU] and O it follows that the spectrum of H is almost surely equal 



to the convex set E + [-B, B] provided B > e 9 + e~ 9 . This is quite different from 
what occurs for the same problem on l 2 (—N,N) in the limit N — > oo, whatever 
boundary conditions are assumed. The point is that for such operators approx- 
imate eigenvalues need not be close to genuine eigenvalues, and the full spectral 
behaviour of the operators can best be seen using pseudospectral theory. 

The determination of the spectrum for smaller values of B is more complicated. 
An application of Theorem [1J] shows that if B < e 9 — e~ 9 then 

Spec(tf) n {z G C : \z\ < r} = 

where r = e 9 — e~ 9 — B > 0. Since E C Spec (if), it follows that the spectrum of 
H has a hole in it. 

We have not been able to determine the precise range of values of g, B for which the 
spectrum contains a hole. A numerical solution seems out of the question because 
the spectrum is heavily affected by extremely infrequent 'regular structure' in the 
potentials. 



7 Classification of the Spectrum 

The classification of the spectrum of non- self-adjoint operators is in a primitive 
state by comparison with that of self-adjoint operators. We start at an abstract 
level. 

If / G Z 2 (Z) satisfies ||/|| = 1, we put 

var(/) = (Q 2 /,/)-(Q/,/) 2 

where Q is the position operator, provided this is finite. If A is a bounded linear 
operator acting on / 2 (Z), we say that A lies in its localized spectrum a\ oc (A) if there 
exists c and a sequence /„ G / 2 (Z) such that ||/ n || = 1 and var(/ n ) < c for all n 
and 

lim \\Afn-Xf n \\ =0. 

n^oo 

A more general definition is possible but not necessary for our purposes. The 
localized spectrum of A includes all eigenvalues whose associated eigenvectors have 
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finite variances, but need not be a closed set. We say that B lies in the translation 
limit class of A if there exists a sequence of unitary translation operators U n such 
that U*AU n converges strongly to B as n — > oo. 

Lemma 14 If B lies in the translation limit class of A then a\ oc (A) D a\ oc (B). In 
particular A G ai oc (A) for every eigenvalue A of B whose corresponding eigenvector 
has finite variance. 

Proof This is the same as that of Lemma 5, with the extra observation that the 
variance is unchanged by translations. 

The following theorem establishes that the localized spectrum is distinct from 
the approximate point spectrum. We say that an operator A acting on l 2 is a 
finite order convolution operator if it is of the form Af = k * f where * denotes 
convolution and k has finite support. 

Theorem 15 If H = A + B where A is a finite order convolution operator and 
B is compact, then its localized spectrum is the set of all its eigenvalues whose 
corresponding eigenvectors have finite variance. If B = this set is empty. 

Proof Let A G (T loc (H) and let f n satisfy ||/ n || = 1, var(/ n ) < c and ||iJ/ n -A/ n || -> 
as n — > oo. If we put a n = (Qf n , f n ) then we have two cases to consider. 

If a n is an unbounded sequence then by passing to a subsequence we may assume 
that a n diverges. If we define U n by 

\U n (p) m 0m— a n 

and put g n = U*f n then (Qg n ,g n ) = and var(# n ) = var(/ n ). Hence g n lies in the 
compact set 

S = {0 G / 2 (Z) : ||0|| = 1 and (Q 2 0, <P) < c} 

for all n. By passing to a convergent subsequence we assume further that g n 
converges to a limit g G S. Now we also have 

lim \\U*HU n g n - Xg n \\ = 

where U*HU n converges strongly to A as n — > oo. Hence Ag = Xg, which con- 
tradicts the fact, proved using Fourier analysis, that the point spectrum of A is 
empty. 

The alternative case is that there exists a such that \a n \ < a for all n. It then 
follows that f n lies in the compact set 

T = {0 G / 2 (Z) : \\(f)\\ = 1 and (Q 2 0, 0) < c + a 2 } 
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for all n. By passing to a convergent subsequence we assume further that f n 
converges to a limit / G T. It is now immediate that Hf = Xf. 

The final special case of the theorem follows from the fact that A has empty point 
spectrum. 

In spite of the above, in the following context the localized spectrum is quite 
different from the point spectrum. We place ourselves in the situation described in 
the first paragraph of Section 6. Using the formula (|J) for the ellipse E, we define 
I z to be the open interior of E + z and O z to be the open exterior of E + z for any 
zeC. 

Theorem 16 If H is defined by (fjj then 

aUH)D (J {i a no p } 

a,(3eM 

almost surely. 

Proof Given a, (3 G M the operator K defined by 

{Kf) n = 



e 9 j n -\ + afn + e 9 fn+i if n > 1 
e-Vn-i + Pfn + e°f n+1 if n < 



lies in the translation limit class of H by Lemma 5, so it is sufficient to prove that 
if A G I a fl Op then A is an eigenvalue of K whose corresponding eigenvector has 
finite variance. 

The eigenvalues and eigenvectors of K may be determined explicitly. For n > 2 
we solve the equation 

e~ 9 u~ l + a + e 9 u = A 

and observe using the Residue Theorem that A G J Q implies that there are two 
solutions Ui, U2 which both satisfy \u{\ < 1. For n < we solve 

e- 9 v- 1 + (3 + e 9 v = A 

and observe similarly that A G Op implies that there is a solution v which satisfies 
\v\ > 1. We now define 



v n if n < 
ci«" + C2M2 if w < 1. 



/n ~ 

By an appropriate choice of ci, c 2 we can ensure that Kf = Xf as required. 
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8 Higher Dimensions 



The results of the last sections can all be extended to higher dimensions, and we 
briefly describe the situation in Z 2 (Z 2 ). We assume that 

H fm,n 6 9 fm—l,n ^ fm+l,n 6 fm,n—l 6 frn,n+l Vin,nfm,n 

where g > 0, h > and V is a random real-valued potential. We assume that V 
has independent values at different points and that they are distributed according 
to a probability measure \i with compact support M C R. If V — then Fourier 
analysis shows that 

Spec(H) = E + F 

where E, F are the ellipses 

E = {e 9+ie + e~ 9 ~ ie : 9 e [0,2tt]} 

and 

F = {e h+ie + e~ h ~~ ie : 6 G [0,2tt]}. 



It is routine to show that Theorem 12 has the analogue 



Theorem 17 With probability one, the spectrum of H contains the set E + F + m 
for every m 6 M . 

Localization of the spectrum is more complicated but the following ideas provide 
some information. Let A, B be two bounded real- valued functions defined on Z 
and define 

(A + B) m,n A m -\- B n . 



Theorem 18 If A + B lies in the translation limit class of V almost surely, then 

o-\ oc (H) D cri oc (Fi) + a] oc (H 2 ) 
almost surely, where Hj act on l 2 iT^) according to 

Hxf n = e' 9 fn~i + e 9 fn+i + A n f n 

and 



H2fn = e h f n -l + e k fn+l + B n f n 
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Proof We first observe that H\® I -\- 1 ® H<2 lies in the translation limit class of 
H almost surely, so 

a loc (H) D a loc {H x ®I + I®H 2 ). 
By considering test functions of the form fx (g) / 2 we also see that 

(7i oc (F 1 ® I + I (g) H 2 ) D a^Hx) + a loc (H 2 ). 



It might seem that the hypothesis of this theorem is rather special. However if 
the support M of \x is an interval then the theorem is of real value. Under these 
assumptions if a, /3 € M then A + B lies in the translation limit class of V almost 
surely, where 

J a/2 if m > 
1 /3/2 otherwise. 



A - R 



The localized spectrum of such operators is considered in Theorem 16. 



9 Conclusions 



We have analyzed the spectral behaviour of a family of randomly generated non- 
self-adjoint matrices by a variety of different methods. Our conclusion is that 
their eigenvalues depend very sensitively on the matrix entries even for quite small 
matrix sizes. The standard proofs of the existence of the eigenvalues depend upon 
the fundamental theorem of algebra, of which there are many proofs, but it is clear 
that the well-known instability of the roots of high degree polynomials leads to 
the large value of the instability index of general non-self-adjoint matrices. This 
problem does not occur for self-adjoint matrices because the variational theorem 
implies that the eigenvalues of such matrices do not change much under small 
perturbations, ||. Our results indicate that nothing of a comparable nature is 
likely to be available in the non-self-adjoint case. 

We have also investigated a family of randomly generated non-self-adjoint bounded 
operators acting on an infinite-dimensional Hilbert space, for which the eigenvec- 
tors almost surely generate a dense linear subspace. In spite of this the eigenvectors 
almost surely do not form a basis and the eigenvalues almost surely generate only a 
small part of the spectrum. We finally proved that the full spectrum of the infinite 
volume non-self-adjoint Anderson model bears little relationship with the infinite 
volume limit of the spectra of the same operator in finite intervals. 
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